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Abstract.  In  the  Dual-Primal  EETI  method,  introduced  by  Farhat  et  al.  [5],  the  domain  is 
decomposed  into  non-overlapping  subdomains,  but  the  degrees  of  freedom  on  crosspoints  remain 
common  to  all  subdomains  adjacent  to  the  crosspoint.  The  continuity  of  the  remaining  degrees  of 
freedom  on  subdomain  interfaces  is  enforced  by  Lagrange  multipliers  and  all  degrees  of  freedom  are 
eliminated.  The  resulting  dual  problem  is  solved  by  preconditioned  conjugate  gradients.  We  give 
an  algebraic  bound  on  the  condition  number,  assuming  only  a  single  inequality  in  discrete  norms, 
and  use  the  algebraic  bound  to  show  that  the  condition  number  is  bounded  by  C(1  -h  \o^{H/h))  for 
both  second  and  fourth  order  elliptic  selfadjoint  problems  discretized  by  conforming  finite  elements, 
as  well  as  for  a  wide  class  of  finite  elements  for  the  Reissner-Mindlin  plate  model. 


1.  Introduction.  This  article  is  concerned  with  convergence  bounds  for  an  iter¬ 
ative  method  for  the  parallel  solution  of  symmetric,  positive  definite  systems  of  linear 
equations  that  arise  from  elliptic  boundary  value  problems  discretized  by  finite  ele¬ 
ments.  The  original  Finite  Element  Tearing  and  Interconnecting  method  (FETI)  was 
proposed  by  Farhat  and  Roux  [9].  The  FETI  method  consists  of  decomposing  the 
domain  into  non-overlapping  subdomains,  enforcing  that  the  corresponding  degrees 
of  freedom  on  subdomain  interfaces  coincide  by  Lagrange  multipliers,  and  eliminat¬ 
ing  all  degrees  of  freedom,  leaving  a  dual  system  for  the  Lagrange  multipliers.  The 
dual  system  was  solved  by  preconditioned  conjugate  gradients  with  a  diagonal  pre¬ 
conditioner.  Evaluation  of  the  dual  operator  involves  the  solution  of  independent 
Neumann  problems  in  all  subdomains,  and  of  a  small  system  of  equations  for  the 
nullspace  component. 

Farhat,  Mandel,  and  Roux  [8]  recognized  that  this  system  for  the  nullspace  com¬ 
ponents  plays  the  role  of  a  coarse  problem  that  facilitates  global  exchange  of  informa¬ 
tion  between  the  subdomains,  causing  the  condition  to  be  bounded  as  the  number  of 
subdomains  increases.  They  also  replaced  the  diagonal  preconditioner  by  a  block  pre¬ 
conditioner  with  the  solution  of  independent  Dirichlet  problems  in  each  subdomain 
and  observed  numerically  that  this  Dirichlet  preconditioner  results  in  a  very  slow 
growth  of  the  condition  number  with  subdomain  size.  Mandel  and  Tezaur  [13]  proved 
that  the  condition  number  grows  at  most  as  log^  (H/h),  where  H  is  subdomain  size 
and  h  is  element  size,  both  in  2D  and  3D.  Tezaur  [19]  proved  that  a  method  by  Park, 
Justino,  and  Felippa  [16]  is  equivalent  to  the  method  of  [9]  with  a  special  choice  of 
the  constraint  matrices,  and  proved  the  log^ (H/h)  bound  for  this  variant.  For  further 
comparison,  see  Rixen  et  al.  [18]. 

Klawonn  and  Widlund  [10]  have  used  preconditioned  conjugate  residuals  to  solve 
a  saddle  problem  keeping  both  the  original  degrees  of  freedom  and  the  Lagrange 
multipliers,  and  obtained  the  asymptotic  bound  log^  (H/h)  using  an  extension  of  the 
theory  of  [13].  The  saddle  point  approach  has  the  advantage  that  approximate  solvers 
can  be  used  for  both  the  Neumann  and  the  Dirichlet  subdomain  problems,  at  the  cost 
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of  of  solving  a  larger  indefinite  problem  instead  of  a  small  positive  definite  problem. 
In  [10],  Klawonn  and  Widlund  have  proposed  new  preconditioners,  proved  uniform 
bounds  for  modifications  of  the  method  for  the  case  of  coefficient  jumps,  which  include 
an  earlier  algorithm  of  Rixen  and  Farhat  [17],  and  provided  further  theoretical  insights. 

The  original  method  of  [9,  8]  does  not  converge  well  for  plate  and  shell  problems, 
and  the  existence  and  the  form  of  the  coarse  space  depend  on  the  singularity  of  the 
subdomain  matrices.  Therefore,  Mandel,  Tezaur,  and  Farhat  [14]  and  Farhat,  Mandel, 
and  Chen  [7]  proposed  to  project  the  Lagrange  multipliers  in  each  iteration  on  an 
auxiliary  space.  An  auxiliary  space  chosen  so  that  the  corresponding  primal  solutions 
are  continuous  on  crosspoints  made  it  possible  to  prove  that  the  condition  number 
does  not  grow  faster  than  log^  (H/h)  for  plate  problems  [14],  and  fast  convergence  was 
observed  for  plate  [7]  as  well  as  shell  problems  [4].  This  method  is  now  called  FETF2. 
For  related  results  for  symmetric  positive  definite  problems,  see  [6,  18]  and  references 
therein. 

The  subject  of  this  paper  is  the  Dual-Primal  FETI  method  (FETFDP),  intro¬ 
duced  by  Farhat  et  al.  [5].  This  method  enforces  the  continuity  of  the  primal  solution 
at  crosspoints  directly  by  the  formulation  of  the  dual  problem:  the  degrees  of  freedom 
on  a  crosspoint  remain  common  to  all  subdomains  sharing  the  crosspoint  and  the 
continuity  of  the  remaining  degrees  of  freedom  on  the  interfaces  in  enforced  by  La¬ 
grange  multipliers.  The  degrees  of  freedom  are  then  eliminated  and  the  resulting  dual 
problem  for  the  Lagrange  multipliers  is  solved  by  preconditioned  conjugate  gradients 
with  a  Dirichlet  preconditioner.  Evaluating  the  dual  operator  involves  the  solution  of 
independent  subdomain  problems  with  nonsingular  matrices  and  of  a  coarse  problem 
based  on  subdomain  corners.  The  advantage  of  this  method  is  a  simpler  formulation 
than  the  methods  of  [14,  7],  there  is  no  need  to  solve  problems  with  singular  matrices, 
and  the  method  was  observed  to  be  significantly  faster  in  practice  for  2D  problems. 
However,  the  design  of  a  good  3D  variant  of  the  method  is  an  open  problem  [5]. 

In  this  paper,  we  prove  that  the  condition  number  of  the  FETI-DP  method  with 
the  Dirichlet  preconditioner  does  not  grow  faster  than  log^{H/h)  for  both  second 
order  and  fourth  order  problems  in  2D.  By  spectral  equivalence,  the  result  for  fourth 
order  problems  extends  to  a  large  class  of  Reissner-Mindlin  elements  for  plate  bending 
as  in  [11,  12,  14].  After  few  initial  definitions,  which  are  substantially  different,  the 
analysis  is  related  to  the  analysis  developed  in  [13,  14].  Just  as  the  formulation  of  the 
present  method  is  simpler,  the  analysis  is  simpler  and  more  elegant  than  in  [13,  14]. 

The  paper  is  organized  as  follows.  The  notation  and  assumptions  are  introduced  in 
Section  2.  In  Section  3,  we  review  the  algorithm  from  [5].  Section  4  gives  the  algebraic 
condition  number  estimates;  these  estimates  apply  to  partitioning  any  symmetric 
positive  definite  system,  not  necessarily  originating  from  partial  differential  equations. 
Finally,  in  Section  5,  we  prove  the  polylogarithmic  condition  number  bounds  for  two 
model  problems  of  the  second  and  fourth  order. 

2.  Domain  partitioning,  notation,  and  assumptions.  We  are  concerned 
with  iterative  solution  of  symmetric,  positive  definite  linear  algebraic  systems.  Parti¬ 
tioning  of  the  system  is  motivated  as  follows.  Let  H  be  a  domain  in  3?^  decomposed 
into  Ns  non-overlapping  subdomains  . . . ,  12'^%  and  each  of  the  subdomains 

be  a  union  of  some  of  the  elements.  Let  u'*  be  the  vector  of  degrees  of  freedom  for 
the  subdomain  corresponding  to  a  conforming  finite  element  discretization  of  a 
second  order  elliptic  problem  or  a  fourth  order  plate  bending  problem  defined  on  H. 
Let  AT®  and  /®  be  the  local  stiffness  matrix  and  the  load  vector  associated  with  the 
subdomain  H®.  We  denote  the  edges  of  the  subdomains  by  F®*  =  912®  n  912*.  Corners 
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are  endpoints  of  edges. 

The  subdomain  vectors  are  partitioned  as 


(2.1) 


u  = 


■u® 

ut 


where  rtf  are  the  values  of  the  degrees  of  freedom  in  the  subdomain  interior,  u®  the 
values  of  the  degrees  of  freedom  at  the  corners  of  the  subdomain,  and  u®  are  the 
remaining  values  of  the  degrees  of  freedom,  i.e.  those  located  on  the  edges  of  the  sub- 
domains  between  the  corners.  The  subdomain  matrices  are  partitioned  accordingly. 


(2.2) 


iL® 

ir 

^^IC 

ri 

iL® 

rc 

Cl 

iL® 

cr 

■^^CC 

We  use  the  block  notation 


U  = 

and  K  =  diag(iL®)  = 

■ 

_  _ 

RNs  _ 

The  block  vectors  Uc  and  Uj.  of  all  internal,  all  corner,  and  all  remainder  degrees 
of  freedom,  respectively,  are  then  defined  similarly. 


'  uj 

ul 

Ui  = 

,  Uf  — 

,  and  Uc  = 

1 

1 _ 

1 

s 

_ 1 

- 1 

_ 1 

Vectors  of  values  of  degrees  of  freedom  on  the  whole  (911®  and  the  corresponding  block 
vectors  will  be  written  as 


r  t”®  1 

5  ^r.c  — 

^r,c 

,  and  also  Vr  c  = 

Vr 

yNs 

^r,c 

Vc 

.  _ 

Let  the  global  to  local  map  be  a  0  —  1  block  matrix 

r  ^ 


L  = 


That  is,  for  a  global  vector  of  degrees  of  freedom  u®,  is  the  vector  of  correspond¬ 
ing  degrees  of  freedom  on  H®.  The  map  L  is  introduced  so  that  we  can  state  the 
problem  to  be  solved, 

(2.3)  LFKLu<^  =  L'^f, 


independently  of  the  solution  algorithm.  Note  that  ImL  is  the  space  of  all  vectors  u 
that  are  continuous  across  the  subdomain  interfaces  and  that  L  is  of  full  column  rank. 

We  assume  that  each  matrix  iL®  is  symmetric  positive  semidefinite  and  that  iL® 
is  positive  definite  on  the  subspace  of  vectors  that  are  zero  at  subdomain  corners, 
{m®|m®  =  0}.  We  also  assume  that  the  global  stiffness  matrix  L^KL  is  positive 
definite,  or,  equivalently,  that  K  is  positive  definite  on  ImL.  These  assumptions  are 
satisfied  in  the  intended  finite  element  applications. 
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3.  Formulation  of  the  algorithm.  In  this  section,  we  review  the  algorithm 
proposed  in  [5]  in  a  form  suitable  for  our  purposes. 

The  degrees  of  freedom  from  both  sides  of  each  edge  T'**  should  coincide, 

(3.1)  ~  u*  Ips*  =  0. 

In  (3.1),  each  pair  of  subdomains  {s,  t}  is  taken  only  once,  with  the  order  (s,  t)  chosen 
arbitrarily.  We  write  the  constraints  (3.1)  as 

BrUr  =  0,  Br  =  [B^. ,  .  .  .  ,  B^^]. 

Note  that  it  follows  immediately  from  the  definition  of  Br  that 

(3.2)  BrB^  =  21 
and,  for  any  edge 

(3.3)  {B^BrUrYlr^^  =  ±«|r-  -  <|r-)- 


Let  Be  be  a  matrix  with  0, 1  entries  implementing  the  global-to- local  map  on 
subdomain  corners.  That  is,  the  equation 


Uc  —  BqU^j  Be 


B 


1 

C 


determines  the  common  values  of  the  degrees  of  freedom  on  subdomain  corners  from 
a  global  vector  uf. 

From  the  construction,  the  space  of  all  vectors  of  degrees  of  freedom  continuous 
across  the  interfaces  can  be  written  as 


ImL  =  {u\BrUr  =  0,Ue  G  ImBe}- 

The  problem  (2.3)  is  reformulated  as  the  equivalent  constrained  minimization  problem 
-u^Ku  —  f  min, 

subject  to  BrUr  =  0  and  Ue  =  BeU^  for  some  m®, 
which  is  in  turn  equivalent  to  finding  the  stationary  point  of  the  Lagrangean 


C{ui,Ur,ul,\)  =  Kv  -  V  f  +  Ur  B^X,  V  =[?;'*],  = 


Ur 

B^eUS 


Eliminating  uf,  u*,  and  from  the  Euler  equations,  we  obtain  a  dual  system  of  the 
form,  cf.  [5,  Eq.  (14)], 


(3.4) 


F\  =  g, 


and  solve  it  using  the  preconditioned  conjugate  gradients  method  with  the  precondi¬ 
tioner 


(3.5)  M  =  BrSrrBj, 
where 

(3.6)  5^^  =  diag(S'^J,  S^rr  =  K^r  - 

For  details  of  the  implementation  and  numerical  results,  see  [5]. 
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4.  Algebraic  bounds.  In  this  section,  we  prove  bounds  on  the  condition  num¬ 
ber  of  the  iterative  method  defined  by  Eqs.  (3.4)  and  (3.5).  Denote  ||m||  =  v  vT'u  and, 
for  a  symmetric  positive  semidefinite  matrix  A,  denote  the  induced  matrix  seminorm 
\u\a  =  Vu'^Au  =  ||A^/^u||.  If  A  is  known  to  be  positive  definite,  we  write  \\u\\a 
instead  or  |m|a)  because  the  seminorm  is  then  known  to  be  a  norm. 

From  the  minimization  property  of  the  Schur  complement,  we  immediately  ob¬ 
tain  the  following  lemma,  which  characterizes  the  bilinear  form  associated  with  the 
preconditioner. 

Lemma  4.1.  It  holds  that  uj SrrUr  =  T[mi{v^ Kv\vr  =  UriUc  =  0}. 

The  next  lemma  gives  a  more  specific  description  of  the  matrix  of  the  dual  equa¬ 
tion  (3.4). 

Lemma  4.2.  It  holds  that  F  =  BrS~^B^ ,  where  the  positive  definite  matrix  S  is 
defined  by 

(4.1)  SUr  =  T[m].{v* Kv\vr  =  Ur,Vc  €  ImBc}- 

Proof  Let 

L{ur,  A)  =  min  L{ui,  Ur,  u®,  A) 

Ui  ,Uc 


Then, 

(4.2)  £{Ur,X)  =  Sur  +  u^bJ X  -  u^hr, 

with  some  hr-  Minimizing  over  Ur,  we  get  Ur  =  S~^{hr  —  BjX).  Substituting  Ur 
into  (4.2)  and  taking  the  variation  over  A  gives  (3.4)  with  F  =  BrS~^Bf .  □ 

We  can  now  characterize  the  norm  induced  by  the  dual  matrix  F. 

Lemma  4.3.  It  holds  that 

(4.3)  = 


Proof.  From  Lemma  4.2, 

||Wr|r 

The  substitution  Wr  =  S^^'^Vr  yields  (4.3).  □ 

Let  y  =  be  the  space  of  Lagrange  multipliers.  In  this  space,  define  the  norm 


and  also  the  dual  norm, 

(4.4)  ||A||v,  =  maxj^^^. 

Since  V  =  ImBr,  substituting  /i  =  Bj-Wr,  we  can  rewrite  (4.4)  as 


(4.5) 


\w. 


^BfX\ 


I  y/  =  max 

\\Bj  Br, 


Wr 
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The  main  result  of  this  section  is  the  following  theorem,  which  gives  a  bound  on 
the  minimal  and  maximal  eigenvalues  of  the  preconditioned  operator  MF. 

Theorem  4.4.  If  there  exists  a  constant  ci  such  that  for  all  Wr, 

\\B'^ BrWrWs^^  <  ClllWrlll, 

then 

>^max{MF)  ^  ^ 

-  4  ■ 


Proof.  The  proof  is  based  on  a  comparison  of  (4.5)  and  (4.3).  Using  Lemma  4.3, 
the  substitution  Vr  =  B^B^Wr,  and  the  property  (3.2),  we  find  that 


X^FX  =  max 

Vr-^O 


>  max 


=  4 


■ic, 


max 


BrWr^O  WBfBrWr 


2  ■ 

S 


Since,  by  definition,  WBJ BrWrW'^g  <  \\Bf  BrWr\\s^^,  we  conclude,  comparing  with  (4.5), 
that 


(4.6) 


A^FA>4||A||^,. 


On  the  other  hand,  from  the  assumption,  we  obtain 


fA  \Ttp\  \wfBf'X\‘^ 

(4.7)  A  UA  =  max — - r-;^ — <  ci  max  ,,  „  ^ —  =  ci| 

Wr.^0  IllUrlll  Br-Wr-^O  \\ Bf  BrWr\\ Sr-r- 

Trivially  from  the  definition  of  M  and  the  norm  in  V,  we  have 

(4.8)  WfiWl,  = 


■'^BfX\ 


lo¬ 


using  (4.6),  (4.7),  and  (4.8)  in  [13,  Lemma  3.1]  completes  the  proof.  □ 

We  now  show  how  to  verify  the  assumption  of  Theorem  4.4  from  inequalities  of  a 
form  that  is  more  usual  in  substructuring  and  easier  to  estimate  for  boundary  value 
problems. 

Denote  by  U®’*  the  operator  that  extends  the  vector  of  values  of  degrees  of  freedom 
on  r®*,  not  including  corners,  by  zero  entries  to  a  vector  of  values  of  degrees  of  freedom 
on  the  whole  (9U®,  and  let  £s  be  the  set  of  all  indices  of  neighbors  U*  of  the  domain 
U®,  with  a  common  edge  T®*.  Denote  by  S'®  the  Schur  complement  on  (dU®  obtained 
by  eliminating  the  interior  degrees  of  freedom  of  U®,  i.e., 

■u®^S®u®  =  min{u®^iL®u®|u®_^  =  u® 

Our  estimate  is  based  on  an  a-priori  bound  of  the  error  of  approximating  a  vector 
of  interface  degrees  of  freedom  that  is  continuous  across  the  corners  by  a  vector  that 
is  continuous  also  across  the  edges.  In  the  applications  in  Sec.  5,  the  approximating 
vector  will  be  chosen  as  the  natural  interpolation  on  the  edges  from  the  corners. 

Theorem  4.5.  Suppose  there  is  a  constant  C2  such  that  for  every  Wr,c,  Wc  € 
Imi?c,  there  exists  Ur  such  that  BrUr  =  0  and,  for  all  s  and  all  t  G  £s, 

(4.9)  |L;®’*(w; -0||,  <  C2|u;;_J|,,  i  =  s,t. 
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Then, 


XmaxiMF) 

XrmuiMF)  - 


where  Ue  is  the  maximum  number  of  the  edges  of  any  subdomain. 

Proof.  Let  Wr  be  given  and  define  Wc  to  be  the  optimal  corner  degrees  of  freedom 
from  the  definition  of  S,  cf.  (4.1).  Then 


Ns 

(4.10)  lk.ll|  =  El<cls- 

Let  uf.  be  as  in  the  assumption  of  the  theorem.  Then  BrUr  =  0,  and,  consequently, 

BJ  BrWr  =  B1^  Br{Wr  —  Mr)- 

Extending  Bf  BrWr  by  zero  values  of  all  corner  degrees  of  freedom,  we  get  using  the 
definition  of  Srr,  cf.,  (3.6),  that 


Ns 

(4.11)  \\BjB,.Wr\\l^=Y,Kc\%, 

S=1 

where 

=  (BJ Br{Wr  —  Mr))*  ,  M®  =  0. 

Using  the  definition  of  E®’*,  we  have 

M®  ,  =  E 

Hence,  from  the  triangle  inequality,  and  then  using  the  property  (3.3)  of  Br  and  the 
triangle  inequality  again,  it  follows  that 

|m®  <  E 

<  E  (|i?*’‘(^r  -  K)\ss  +  -  m‘)|s0  • 

t^E  g 

Squaring  the  the  last  inequality,  using  the  inequality  (a  +  b)^  <  2(a^  +  &^),  and  the 
a-priori  bound  (4.9)  yields 

^  2c2  E  (l^r,cll=  +  l^r.cll*)  • 

t^Es 

By  the  summation  over  the  subdomains  and  using  (4.11)  and  (4.10),  we  can  conclude 
that 


\BT  BrWr 


Ns 

El 

S=1 


<  4c2n, 


Ns 

XI 

s=l 


=  4c2ne||M;r 


Is 


It  remains  to  use  Theorem  4.4.  □ 
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5.  Applications.  In  this  section,  we  verify  the  assumption  of  Theorem  4.5  for 
two  model  problems.  The  Sobolev  seminorm,  denoted  by  is  the  U‘{X)  norm 

of  the  generalized  derivatives  of  of  order  to  >  0  of  the  function  u.  The  Sobolev 
norm  is  then  defined  by  ||M||m,p,x  =  ||  [|m|o,p,X)  ■  •  ■ ,  \u\m,p,x\  \\ip-  Sobolev  norms  for 
noninteger  to  are  defined  by  interpolation.  Cf.,  e.g.,  [2,  15]  for  details  and  references. 
In  particular,  jujopx  =  ||m||opX  =  ||u||lp(x))  and,  on  the  boundary  T  of  a  domain  in 


\u{x)  - 


dx  dy. 


r  JT 


5.1.  A  second  order  elliptic  problem.  Consider  the  boundary  value  problem 


(5.1) 

where 


Au  =  g  in  n. 


M  =  0  on  did, 


with  a{x)  a  measurable  function  such  that  0  <  oq  <  a{x)  <  ai  a.e.  in  il. 

We  assume  that  model  problem  (5.1)  is  discretized  using  conforming  PI  or  Q1 
elements  and  denote  by  the  corresponding  finite  element  space  that  satisfies 

the  usual  regularity  and  inverse  properties  [2] .  Let  h  denote  the  characteristic  element 
size.  Assume  that  all  functions  in  vanish  on  the  boundary  of  il.  Assume 

that  each  subdomain  is  the  union  of  some  of  the  elements  and  denote  the  space  of 
restrictions  of  functions  from  to  subdomains  by  ^^^(n®).  For  every  vector 

of  degrees  of  freedom  m®,  denote  by  /piM®  the  corresponding  finite  element  function. 
The  trace  of  this  function  is  determined  by  degrees  of  freedom  on  the  boundary  only 
and  we  often  abuse  the  notation  to  define  the  trace  from  the  boundary  values  only. 

For  simplicity,  we  assume  that  the  subdomains  Hi,  i  =  l,...,Ns  form  a  quasi¬ 
regular  triangulation  of  the  domain  il.  Denote  the  characteristic  size  of  the  subdo¬ 
mains  by  H.  Finally,  let  C  denote  a  generic  constant  independent  of  h  and  H . 

Theorem  5.1.  For  the  model  second  order  problem,  it  holds  that 


>^max{MF) 

XrmniMF) 


<C[l  +  logj 


where  the  constant  does  not  depend  on  FI  and  h. 

Proof.  We  only  need  to  verify  the  assumptions  of  Theorem  4.5.  For  a  given  w^^c-, 
define  Ur^c  by  linear  interpolation  from  Wc',  that  is,  Ipitt®  ^  is  linear  on  all  edges  of  fl®, 
and  Uc  =  Wc-  Writing  Wr,c  =  {wr,c  —  Ur^c)  +  Ur^c,  we  get  from  [1,  Lemma  3.5]  that 


1/piF;®’‘(u;®-<)1 


i,2,an=  ^  C' 


1  -f  log 


Using  the  uniform  equivalence  of  the  seminorms  \Ipivff\x ^2,dQ,=  ~  br,cls“)  cL,  e.g., 
[1],  and  the  uniform  equivalence  of  the  seminorms  juj  i  2  90=  ~  1^1 5  2  an*  ^or  functions 
with  support  on  the  edge  F®‘,  which  follows  from  the  fact  that  the  subdomains  are 
shape  regular,  we  get  (4.9)  with  C2  =  C  (l  -I-  log  .  □ 

Remark  5.2.  The  assumptions  that  the  subdomains  have  straight  edges  can  be 
relaxed  to  accommodate  the  case  of  only  shape  regular  subdomains  with  edges  that  are 
not  straight  in  a  standard  way  by  mapping  from  one  or  more  reference  subdomains. 


5.2.  A  fourth  order  problem.  Consider  a  biharmonic  boundary  value  problem 
in  a  variational  form:  Find  u  G  such  that 

(5.2)  a{u,v)  =  f{v),  yvGHg{n), 
where 

a{u,v)  =  /  diiudiiv  +  di2udi2V  +  d22ud22V,  yu,v  G 

Jn 

Let  the  model  problem  (5.2)  be  discretized  by  reduced  HCT  elements  [2].  We  use 
the  same  assumptions  on  the  decomposition  as  in  Section  5.1;  in  particular,  the  sub- 
domain  edges  are  straight.  Let  V'^^^(n)  be  the  finite  element  space  of  HCT  elements 
satisfying  the  usual  regularity  and  inverse  properties  and  the  essential  boundary  con¬ 
ditions.  Note  that  C  (7^(12)  n  iLg(r2).  On  each  element,  a  function  v  in 

is  determined  by  the  values  v{ai)  and  the  values  of  its  derivatives 
j  =  1,2,  at  the  vertices  of  the  element.  Denote  by  t,  n,  dt,  and  i9„  the  tangential 
and  normal  directions,  and  the  tangential  and  normal  derivative,  respectively.  The 
traces  of  functions  from  V)f^^^(O)  on  50  are  pairs  of  functions  {u,  Vu)  such  that  u  is 
piecewise  cubic,  n  •  Vu  is  piecewise  linear,  and  u  and  Vu  are  consistent,  dtu  =  t  -Vv. 
We  will  abuse  notation  and  write  the  trace  functions  simply  as  (m,  Vu)  or  just  u.  The 
space  of  HCT  trace  functions  is  denoted  by  Denote  the  finite  element 

interpolation  operator  by  Ihct-  As  in  Section  5.1,  we  abuse  the  notation  by  defining 
the  trace  of  the  interpolation  from  the  boundary  degrees  of  freedom  and  write,  e.g., 
lHCTu’^r,c  y  We  adopt  the  convention  that  all  functions  are  understood 

extended  by  zero  outside  of  their  stated  domain;  e.g.,  u|r  is  zero  outside  of  F. 

The  following  lemma  gives  a  bound  on  the  interpolation  from  subdomain  corners 
to  edges. 

Lemma  5.3.  For  every  w  G  V'^'"^(5r2®),  define  u  G  {dVl^)  by  HCT  in¬ 

terpolation  from  the  corners  of  to  the  edges;  that  is,  u  =  v  and  Vu  =  Vw  at  the 
corners  of  ,  and,  on  each  edge  F®*,  u  is  a  cubic  polynomial  and  n  ■  Vu  is  linear. 
Then 

(5.3)  |(^(^-^)lr=*)|i,2,an=  ^  ^ 

Proof.  Denote  v  =  w  —  u  G  {dPt).  Then  r;  =  0  and  Vu  =  0  on  the  corners 

of  n®,  so  from  [12,  Lemma  4.1],  it  follows  that 

(5.4)  |(VHr-)|5,2,aa=  ^  ^  (l  +  log  f )  l|Vu||2,^,r- 

Since  Vwlrs*  is  from  a  space  of  dimension  6  and  all  norms  on  a  finite  dimensional 
space  are  equivalent,  we  have  in  the  case  when  the  length  of  F®*  is  one  that 

^  li^'^ll i,2,r'’‘  ~  II Vujjo.oo.r''*. 

Since  |Vw|  i  2,r''‘  ||  Vwllo^oc.r"*  are  invariant  to  stretching  the  edge,  we  get 

(5.5)  |Vu|i,2.r- <  ^^l|Vw||o,oo,r- 

in  the  general  case  by  scaling.  We  will  show  at  the  end  of  the  proof  that 

(5.6)  llVullo.oo.r-  <5||Vu;||o.oo.r-. 
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From  the  discrete  Sobolev  inequality  as  generalized  by  [12,  Lemma  4.2],  we  have 

(5.7)  <  C*  (^1  +  log  (^|Vw||_2.r=‘  +  ;^l|V^«llo.2.r=‘^ 
Combining  (5.4),  (5.6),  and  (5.7),  we  obtain 

(5.8)  |Vf||_2.on=  ^  C*  (l  +  logf-) 


First  consider  the  case  when  n  d^l  =  0.  Then  (5.3)  is  invariant  to  adding  a  linear 
function  to  w  because  v  is  the  error  of  linear  interpolation  of  w  and  because  this  adds 
only  a  constant  to  Vw.  So,  without  loss  of  generality,  let  \7w  =  0.  Then  by  the 
Poincare-Friedrichs  inequality 


(5.9) 


1 

H 


l!y|lo.2,3n 


^<C\y\ 


2 

5,2,90=’ 


proved  first  on  a  reference  domain  and  then  scaled  [3,  20].  Now  (5.9)  and  (5.8) 
give  (5.3).  In  the  case  when  C  d^l  yf  0,  there  are  some  essential  boundary  condi¬ 
tions  on  Since  (5.3)  has  already  been  proved  in  the  absence  of  essential  boundary 
conditions  on  d^l^,  it  is  sufficient  to  restrict  (5.3)  onto  the  subspace  defined  by  the 
boundary  conditions,  noting  that  u  satisfies  the  boundary  conditions  as  well. 

It  remains  to  prove  (5.6).  The  inequality  is  trivial  for  n  ■  Vm  because  the  normal 
derivative  is  interpolated  linearly  between  the  corners.  Let  ul  be  the  linear  function 
on  the  edge  F^*  defined  by  the  values  of  u  on  the  corners.  Then,  using  the  triangle 
inequality. 


(5.10)  \\dtu\\  o.oo.ro*  <  \\dt{u-UL)\\  o,oo,r^*  +  \\dtUL\\  o,oo,r'^*  ■ 

By  a  simple  computation,  we  see  that  the  Hermite  basis  function  (p,  (()(0)  =  0,  <()'(0)  = 
1,  0(1)  =  0'(1)  =  0,  0  a  polynomial  of  order  3,  attains  the  maximum  of  j0'j  at  0. 
Mapping  the  interval  (0, 1)  on  the  edge  F®*  and  noting  that  u  —  ul  is  zero  at  the 
endpoints  xi  and  X2  of  F®‘,  we  get  using  the  triangle  inequality  that 

\\dt{u-UL)\\  o,oo,r=*  <  \dt{u  -  ul){xi)\  +  \dt{u  -  ul){x2)\ 

(5.11)  <  2j|9tu;llo,oo,r=*  +  2l|9tMLllo,oo,r=*, 

because  dtu{xi)  =  dtw{xi),  i  =  1,2.  From  the  mean  value  theorem  and  the  fact  that 
UL{xi)  =  w{xi),  i  =  1,2,  it  follows  that 

(5.12)  II II o,oo,r=*  <  II 9(^11  o,oo,r=*, 
hence  (5.11)  gives 


(5.13)  \\dt{u  -  ML)||o.oo.r=*  <  4||i9tw||o,oo.r=‘- 

Now  (5.10),  (5.12),  and  (5.13)  give  (5.6).  □ 

Theorem  5.4.  For  the  fourth  order  model  problem, 


XmaAMF) 

X^,n{MF) 


<  C  (  1  -f  log  ^ 


where  the  constant  does  not  depend  on  FI  and  h. 
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Proof.  The  proof  follows  immediately  from  Theorem  4.5,  with  (4.9)  being  a  conse¬ 
quence  of  Lemma  5.3  and  the  uniform  equivalence  ofseminorms  \V  Ihctw^.  2  ~ 

KJs=,cf.,  [12].  □ 

Remark  5.5.  The  result  extends,  by  spectral  equivalence,  to  DKT  elements  and 
a  certain  class  of  non-locking  elements  for  the  Reissner-Mindlin  plate  model  as  in 
[12].  The  result  also  extends  to  the  case  when  the  subdomain  edges  are  not  straight  by 
considering  the  subdomains  to  be  images  of  a  reference  domain. 
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